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PENETRATION  MECHANICS  OF  TEXTILE  STRUCTURES 
I.  BALLISTICS  OF  TRANSVERSELY  IMPACTED  FIBERS* 

Introduction 

Although  impact  of  single  fibers  or  fiber  assem¬ 
blies  is  an  important  subject  in  its  own  right,  being 
relevant  to  climbing  ropes,  aircraft  carrier  arrest 
cables,  high-speed  weaving,  etc.,  the  principal  develop¬ 
ments  in  this  area  have  been  made  by  workers  whose 
major  interests  have  been  in  the  impact  resistance  of 
woven  or  non-woven  textile  structures.  The  most  notable 
of  these  structures  have  been  the  lightweight  armor  vests 
used  by  police  and  military  personnel,  but  among  other 
important  applications  can  be  listed  aircraft  engine 
containment  shrouds,  fl3k  blankets,  and  vehicle  seat 
belts.  Ballistic  nylon  has  been  used  successfully  for 
these  vests  since  the  Second  World  War,  although 
current  developments  have  emphasized  the  du  Pont  aramid 
fiber  marketed  as  Kevlar**.  Although,  as  will  be  shown 
below,  excellent  sing.le-fiber  ballistic  response  does 
not  necessarily  guarantee  a  superior  vest,  any  under¬ 
standing  of  textile  structure  ballistics  must  be  pre- 
*  From  Reference  1  (see  Page  109).  Used  by  permission 
of  Textile  Research  Journal 
**Trademark  of  E.I.  du  Pont  de  Nemours  &  Co.,  Inc. 
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ceeded  by  an  understanding  of  single-fiber  response. 

A  strong  motive  for  discussing  fibers  is  that 
single  fiber  tests  are  often  used  as  screening  tests 
for  ballistic  protection  materials.  As  an  example, 
one  often  encounters  tabulations  of  "transverse 
critical  velocity",  that  ballistic  velocity  at  which 
a  transversely  impacted  yarn  experiences  nearly 
instantaneous  failure.  Typical  data  is  shown  below. 

Transverse  critical  velocities  of  textile 
fibers. I2] * 


Vcr,  m/sec 


Nylon 

616 

Polyester 

472 

Nomex 

442 

Fiberglass 

274 

Kevlar  29 

570 

Such  tests  are  often  indicative  of  relative 
ballistic  resistance,  but  perfect  correlations  cannot 
be  guaranteed.  In  the  above  tabulation  Kevlar  29  proves 
to  be  the  best  ballistic  material  when  put  into  a  panel, 
in  spite  of  its  having  a  lower  transverse  critical 

*  Numbers  in  brackets  refer  to  references  listed  on 
pages  109-U0. 
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velocity  than  nylon. 


Longitudinal  Wave  Propagation 

Wave  propagation  phenomena  in  fibers  and  thin 
rods  are  considerably  less  complicated  than  in  a 
general  medium,  since  the  possibility  of  unrestrained 
transverse  contraction  in  fibers  eliminates  (to  a 
good  approximation)  the  simultaneous  propagation 
of  independent  dilatational  and  distortional  waves 
which  are  present  in  general.  The  equation  of  motion 
for  fibers  or  rods  is  simply  [3] : 


where  u  is  the  longitudinal  particle  displacement,  p 
is  the  material  density,  E  is  the  longitudinal  Young's 
modulus,  and  x  and  t  are  the  space  and  time  coordinates. 
This  is  the  well-known  wave  equation,  whose  solution 
represents  a  disturbance  traveling  at  a  velocity 

C*=.  V&Tp  (2) 
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Conventional  textile  units  employing  stiffness  per  unit 
linear  density  are  very  convenient  in  wave  propagation 
analyses,  since  the  factor  p  is  included  implicitly 
in  the  modulus.  For  modulus  expressed  in  grams  per 
denier  and  wavespeed  in  meters  per  second.  Equation  3 
becomes: 

c  =  Vki.  <3> 

where  k  =  88,260  is  the  necessary  unit s-con version 
factor.  In  these  equations,  as  well  as  those  to 
follow,  the  modulus  is  taken  to  be  the  "dynamic” 
stiffness  relevant  to  the  high  strain  rates  corres¬ 
ponding  to  wave  propagation  tests.  The  development  of 
such  dynamic  constitutive  relations  from  experimental 
fiber-impact  data  has  been  described  elsewhere  [4,5]. 

Consider  a  fiber  fixed  at  one  end  whose  free  end 
is  suddenly  subjected  to  a  constant  outward  velocity 
V  in  the  longitudinal  (fiber)  direction.  After  a  time 
it  ,the  strain  wave iwill  have  propagated  into  the  fiber 
a  distance  ct,  while  the  free  end  will  have  displaced 
outward  an  amount  Vt-.  The  strain  resulting  from  the 
impact  is  then  the  displacement  Vt  divided  by  the 
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affected  length  ct: 


G  =■ 


vt  v 


l/ke. 


(4) 


The  corresponding  stress  is 


er-Es-  VVETIZ 


(5) 


te.  * 


The  above  relations  have  assumed  a  linear 
elastic  material  whose  stiffness  E  is  independent  of 
the  strain.  In  this  case  the  wavefront  will  propagate 
as  a  sharp  discontinuity  (a  shock  wave)  at  which  the 
strain  rises  instantaneously  from  zero  to  the  value 
given  by  Equation  4.  Many  ballistic  fibers  are 
nonlinear,  however,  and  the  effect  of  material  non¬ 
linearity  leads  to  some  complication  of  the  above 
description.  A  nonlinear  fiber  can  be  characterized 
as  having  a  strain-dependent  modulus  E  =  E(  e.  ) ,  so 
that  Equation  3  becomes: 


as 

I  38 


:•  M 


c-  cC6^  «= 


■\/kEl  C£>> 


(6) 
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where  £q  is  the  ultimate  value  of  strain  generated 


by  the  impact.  Since  the  particle  velocity  must  match 


the  imposed  velocity,  we  have 


V-  J 


The  strain  Cn  developed  by  longitudinal  impact  is 


found  by  solving  Equation  8,  perhaps  numerically. 


Transverse  Impact  of  Fibers 


As  the  transverse  impact  of  fibers  seems  intuitive¬ 


ly  germane  to  impact  of  woven  textile  panels,  the 


technical  community  interested  in  lightweight  ballistic 


protection  has  devoted  intensive  effort  to  this  problem 


since  World  War  II.  Following  the  pioneering  works  of 


Taylor  [6]  and  von  Karman  [7]  during  the  war,  valuable 


contributions  have  been  made  by  Peterson  et  al.  [8], 


Shultz  et  al.  [9] ,  Wilde  et  al.  [5] ,  among  others,  but 


by  far  the  most  prolific  of  these  efforts  has  been 


that  of  Jack  C.  Smith  and  his  colleages  at  the  National 


Bureau  of  Standards.  Reference  [10]  provides  a  review 


of  most  of  this  work,  which  contains  a  wealth  of 


&/■  >/*  '  -■ 


agga 
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experimental  and  theoretical  contributions  ranging  over 
a  period  of  approximately  ten  years  in  the  fifties 
and  sixties. 

v 


Figure  1.  Wave  Propagation  in  a  transversely  impacted  fiber. 

The  rate-independent  theory  of  transverse  fiber 
impact  as  developed  by  Smith  can  be  stated  with  refer¬ 
ence  to  Figure  1.  This  illustrates  a  fiber,  originally 
straight  in  the  horizontal  direction,  which  has  been 
impacted  by  a  projectile  traveling  vertically  upward. 
Upon  impact,  longitudinal  waves  of  the  type  described 
in  the  previous  section  are  propagated  outward  from  the 
point  of  impact.  Behind  these  waves  material  flows 
inward  toward  the  point  of  impact  at  a  constant 
velocity  w,  strain  €  and  stress  C^. 


In 


addition  to  the  longitudinal  waves,  transverse  "kink" 
waves  are  also  propagated  outward  from  the  impact  point. 
At  the  transverse  wavefront  the  .inward  material  flow 
ceases  abruptly  and  is  replaced  by  a  transverse 
particle  velocity  equal  in  magnitude  and  direction  to 
that  of  the  projectile.  The  strain  and  tension  are  un¬ 
changed  across  the  transverse  wavefront,  but  both  the 
longitudinal  and  transverse  particle  velocities  ex¬ 
perience  discontinuities  there;  ir.  this  sense  the  trans¬ 
verse  wave  is  a  geometrical  shock.  The  apprently  un¬ 
balanced  tensions  on  either  side  of  the  transverse 

wavefront  are  compensated  by  the  change  in  particle 
momentum  as  the  wave  propagates.  Behind  the  trans¬ 
verse  wavefront  all  particle  velocities  are  equal  in 
magnitude  and  direction  to  the  projectile  velocity, 
and  the  fiber  configuration  is  a  straight  line  at  a 
constant  inclination  0  from  the  longitudinal  direction. 

The  inward  particle  velocity  is  found,  as  in  the 
longitudinal  case,  as 


The  final  strain  c Q  is  unknown  as  yet,  but  E(  e  ) 
is  known  as  the  slope  of  the  dynamic  stress-strain 
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curve.  The  outward  velocity  U  of  the  transverse  kink 
wave,  measured  relative  to  a  Lagrangian  frame  attached 
to  and  extending  with  the  fiber,  is: 


(10) 


To  a  fixed  observer  the  transverse  wave  appears  to 
propagate  in  a  "laboratory"  frame  of  reference  at 

(1=  (i+O  (J  -to  (id 

Finally,  the  above  variables  are  related  to  the  impact 
velocity  V  through  the  relation: 


in  yo+€*f  uz-  u; 


(12) 


Equations  9-12  constitute  four  relations  between,  V, 
w,  €:  q,  cr^,  u'  anc^  The  material  dynamic  stress 
strain  curve  relates  cj*0  and  £  g'  so  t^iat  once  one 
of  the  parameters  (say  V)  is  specified,  the  other  four 
independent  parameters  (w,  e  q,  U,  U)  can  be  found. 
For  nonlinear  stress-strain  curves,  numerical  solu¬ 
tions  will  likely  be  more  convenient. 


Certain  limitations  to  the  Smith  analysis 
described  above  must  be  mentioned.  First,  it  is 


rate-independent.  Most  polymeric  fibers  exhibit 
strong  rate  dependencies,  and  these  effects  are  beyond 
the  capacity  of  this  analysis  to  describe.  Perhaps 
a  more  severe  limitation  is  that  the  Smith  analysis  is 
not  applicable  to  late-time  effects  in  the  wave  propa¬ 
gation  process.  In  real  situations  the  outgoing 
longitudinal  wave  soon  collides  with  an  obstacle:  a 
clamp,  in  the  case  of  single-fiber  tests,  or  a  fiber 
crossover,  in  the  case  of  impact  in  woven  textile 
panels.  Upon  such  a  collision  a  reflected  wave  is 
propagated  from  the  collision  point  in  the  direction 
opposite  that  of  the  original  wave.  This  reflected 
wave  in  turn  soon  collides  with  the  outward -traveling 
transverse  wave,  and  this  collision  generates  another 
two  waves  which  travel  away  from  the  collision  point. 
These  waves  in  turn  eventually  collide  with  the  clamps, 
or  the  projectile,  or  other  waves.  The  result  of 
4  hese  wave  reflections  and  interactions  is  a  situation 
which  becomes  intractable  by  closed-form  mathematical 
methods,  and  this  late-time  intractability  is  a 
principal  reason  for  the  development  of  numerical 
computer  solutions. 

Use  of  the  Rate-Independent  Theory  in  Preliminary  Design 

In  spite  of  the  limitations  of  the  Smith  theory 
outlined  above,  the  rate-independent  analysis  provides 
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a  highly  useful  means  of  assessing  approximate  relations 
between  fiber  material  properties  and  ballistic 


response.  These  relations  are  of  considerable  value 
in  performing  preliminary  design  steps  in  development 
of  textile  ballistic-protection  devices. 

Assuming  the  material  to  be  linear  in  stress- 
strain  response  (E  =  constant) ,  the  Smith  analysis  can 
be  cast  in  the  simple  form: 


V  «  -)J e.0lcE  ]  (13) 


which  provides  a  relation  for  the  strain  €1  Q  developed 
by  impact  at  a  velocity  V  in  terms  of  the  fiber  modulus. 
The  relation  can  be  solved  numerically  if  one  wishes  to 
compute  e:  0  for  a  given  V,  or  it  can  be  used  directly 
to  plot  e  0  versus  V  for  the  purpose  of  developing 
design  curves  (see  Fig.  2) .  Once  €  q  is  known,  then 
o~  q»  U,  U,  and  w  can  be  found  from  either  the  stress- 
strain  relation  or  Equations  9-12.  Figure  3  shows 
such  a  plot  of  tension  ,  (3~q  versus  V,  with  modulus  as  a 
parameter. 
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Figure  2.  Predicted  impact  strain  for  linear  rate-independent  fibers. 
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Figure  3.  Predicted  impact  tension  for  linear  rate-independent  fibers. 
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Since  the  above  curves  rise  monotonically  with 
velocity,  one  can  observe  the  influence  of  modulus  more 
easily  by  plotting  ballistic  response  at  a  constant 

velocity,  and  Figure  shows  such  a  plot  at  V  =  400 
m/sec.  Here  are  plotted  the  strain  and  tension  from 
the  above  methods,  along  with  the  strain  energy 
^  =  j  developed  behind  the  wave  and 

the  rate  of  energy  absorption  /c  of  the  fiber. 

(The  term,  S*c,is  shown  in  mixed  units,  but  it  could 
be  converted  to  joules/sec  once  the  density  and  denier 
of  the  fiber  are  specified.)  The  rate  of  energy 
absorption  at  the  wavefront  must  equal  the  rate  at 
which  the  fiber  extracts  kinetic  energy  from  the 
projectile,  and  it  is  a  reasonable  measure  of  ballistic 
efficiency.  Note  that  this  energy  absorption  rate 
rises  monotonically  with  fiber  modulus,  although  with 
less  dramatic  improvements  after  approximately  500  g/den. 
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BALLISTIC  RESPONSE,  relative  units 


Figure  4*  Effect  of  fiber  stiffness  on  ballistic  response;  10  =  10 
g/den  for  tension,  10$  for  strain,  0.03  g/den  for  strain 
energy,  and  900  g/den  sec  for  energy  absorption  rate. 


Of  course,  one  cannot  improve  ballistic  efficiency 
indefinitely  by  continuing  to  seek  stiffer  fibers.  In 
general,  increases  in  stiffeness  are  accompanied  by  de¬ 
creases  in  breaking  strain,  and  a  point  may  well  be 
reached  where  this  reduced  ductility  overshadows  the 
beneficial  reduction  in  impact-generated  strain  shown  in 
Figure  2.  This  effect  may  be  quantified  by  means  of 
Equation  13,  where  one  may  calculate  the  critical  trans¬ 
verse  velocity  by  determining  the  velocity  which  just 
generates  the  dynamic  breaking  strain  on  impact.  If  one 
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knows  the  variation  of  breaking  strain  with  stiffness, 
these  calculations  may  be  used  to  select  approximately 
an  optimum  fiber  stiffness  for  ballistic  efficiency. 
This  process  is  carried  through  for  illustrative  pur¬ 
poses  in  Figure  5. 


Modulus,  gpd 

Figure  5«  Prediction  of  optimum  stiffness  for  nylon  fibers. 


The  dashed  line  in  this  Figure  is  the  relation 
between  dynamic  stiffness  and  breaking  strain  as 
determined  from  fiber-impact  tests  on  a  series  of 
nylon  yarns  that  had  been  subjected  to  various 


drawing  treatments  by  the  manufacturer  [5] .  The  solid 
line  is  the  the  calculated  transverse  critical 
velocity,  considering  the  effect  of  stiffness  on  both 
impact-induced  strain  and  on  breaking  strain.  An 
optimum  is  observed  near  60  g/den,  which  is  in  fair 
agreement  with  experimental  observation.  All  this  is 
a  quantification  of  the  often-quoted  guideline  in  armor 
design  that  one  seeks  the  highest  possible  modulus  in 
order  to  spread  the  impact  over  a  wide  area  via  in¬ 
creased  wavespeed,  but  that  the  process  must  not  be 
carried  so  far  as  to  induce  excessive  brittleness.  In 
hard  armor  this  reasoning  has  led  to  the  use  of 
ceramic  faceplates  to  give  high  wavespeed,  backed  by 
a  fiberglass  laminate  to  provide  the  needed  toughness. 

Selection  of  a  Failure  Criterion 

The  use  of  simple  ultimate  breaking  strain  as  a 
failure  criterion  in  the  above  example  is  overly 
simplistic,  since  it  does  not  incorporate  the  strong 
temperature  and  rate  dependencies  that  are  known  to 
exist  in  polymeric  material.  A  versatile  fracture  model 
which  does  incorporate  these  dependencies  and  is  still 
computationally  convenient  is  that  due  to  Zhurkov  {,11]  , 
which  states  that  the  lifetime  t  of  a  solid  subjected 
to  a  uniaxial  stress  c5  is  of  the  form 


{  CU’-y'VVET}  (14) 

where  ^  is  a  pre-exponential  factor  with  units  of 
time,  u*  is  an  apparent  activation  energy  for  the 
fracture  process,  is  a  factor  with  units  of 

volume,  R  is  the  gas  constant  (8.314  J/mole  °K) ,  and 
T  is  the  absolute  temperature.  For  constant  tempera¬ 
ture,  Equation  14  reduces  to 


where 


r  “  CL  C  '&0-V> 

e*r  CaVem"} 

p>  =  t*  /  izt 


When  stress  and  temperature  vary  during  the  loading 
process,  one  can  assume  linear  superposition  and  write 
the  Zhurkov  criterion  in  the  form: 


J  e^Cu^-A-VlZT^  =1  an 


In  a  constant-stress-rate  experiment  at  constant 
temperature,  for  instance, 
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(18) 


o-(tVi)  Ct 

/  A±/  {*.  e*p  C-  p.Ck')  [  =  A 

c> 

[wj?  =1 

t  -  ta.  C  l  +  ^c  ^  ^  c 

To  illustrate  the  order  of  rate  dependency  provided  by 

Zhurkov's  model,  Figure  6  shows  a  plot  of  Equation  18 

for  the  case  of  drawn  nylon  fibers.  In  this  figure 

-  2.20  x  1019  sec,  and  =5.13  (g/den)  1 ; 

these  are  the  values  obtained  by  Zhurkov  [11]  by 

fitting  Equation  15  to  creep-rupture  data.  Such 

a  plot  can  be  used  to  depict  the  time-to-break  for  a 
fiber,  and  the  tenacity-at-break,  as  a  function  of 
the  loading  rate. 
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Figure  6.  Variation  of  breaking  tenacity  with  loading  rate  -  Khurkov 
model. 

As  a  more  direct  example  of  the  utilization  of 
Zhurkov's  model  in  fiber  ballistics,  the  stress  pre¬ 
dicted  by  the  Smith  theory  for  a  given  impact  velocity 
and  fiber  modulus  can  be  used  in  Equation  14  to  pre¬ 
dict  the  time  after  impact  at  which  the  fiber  will 
rupture.  This  analysis  predicts  that  there  is  no 
unique  critical  transverse  velocity,  but  rather  a 
range  of  velocities  over  which  the  fiber  will  fail  in 
experimentally  observable  times.  Figure  7  shows  the 
predicted  results  for  drawn  nylon  fibers,  using  an 
assumed  dynamic  modulus  of  80  g/den  with  the  same 
values  of  cL  and  used  in  Figure  6.  This  figura 


Time  to  Break,  microseconds 
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shows  that  at  velocities  above  approximately  775  m/sec, 
rupture  occurs  in  less  than  fifty  microseconds  and 
would  be  counted  in  most  high-speed  photographic 
records  as  having  occurred  instantaneously  upon  impact. 
The  times -to -break  become  exponentially  longer  at 
lower  velocities,  and  failure  will  occur  at  the  clamp 
due  to  wave  reflection  at  times  dependent  on  the 
wavespeed  and  fiber  length.  This  variation  in  what 
may  be  termed  critical  velocity  for  impact  may  make  up 
a  large  part  of  the  scatter  observed  experimentally  in 
determining  critical  transverse  velocities. 


Figure  7.  Variation  in  transverse  critical  velocity  due  to  fracture  rate 
effects. 


An  important  advantage  to  Zhurkov's  model  is  that 


it  is  derivable  in  terms  of  basic  reaction-rate 
fracture  analysis.  As  such,  it  provides  a  means 
whereby  the  materials  scientist  can  predict  materials 
and  processing  modifications  so  as  to  manipulate  the 
fracture  parameters  and  improve  ballistic  performance. 

A  recent  review  [12]  describes  the  basic  implications 
of  reaction  rate  models  such  as  Zhurkov's,  as  well  as 
their  limitations  and  experimental  corroboration. 

The  development  of  these  models  is  somewhat  controver¬ 
sial,  with  several  quite  divergent  approaches  having 
strong  advocates.  Zhurkov's  model  in  particular  is 
often  criticized  as  being  simplistic,  but  is  convenient 
for  use  in  impact  by  virtue  of  its  computational  con¬ 
venience  and  its  ability  to  model  a  wide  range  of 
materials  behavior,  if  only  phenomenologically.  Finally, 
it  should  be  cautioned  that  the  experiments  Zhurkov  used 
in  corroborating  his  model  were  no  faster  than  the  milli¬ 
second  time  scale,  some  three  orders  of  magnitude 
slower  than  ballistic  impact  fractures.  Such  an  extra¬ 
polation  is  clearly  dangerous  and  should  be  verified  by 
additional  experimentation.  The  plot  given  in  Figure 
7  is  in  reasonable  but  not  excellent  agreement  with 

experimental  data  given  in  Smith's  papers,  indicating 
that  the  approach  is  promising  but  needing  of  further 
corroboration . 
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II.  NUMERICAL  ANALYSIS  OF  IMPACT  ON  WOVEN  PANELS 


The  fabric  of  dimension  L  by  L  shown  in  Figure  8 
is  modeled  as  a  network  of  interconnected  fiber 
elements  impacted  at  zero  obliquity  by  a  rigid  missile 
of  mass.M  with  an  initial  striking  velocity  V  .  The 
network  model  rather  than  a  continuous  membrane  is 
employed  here,  for  it  is  not  only  more  consistent 
with  the  discrete  fabric  structure  but  also  leads  to 
better  agreement  with  the  physical  deformation  con¬ 
figuration  (pyramidal  cone  rather  than  hemisphere)  as 
observed  by  Wilde  [13]  in  high  speed  photographs. 

The  elastic  continuum  supporting  the  fabric  deformation 
is  generally  very  flexible  in  transverse  direction: 
bending  effects  are  neglected.  The  constituent  fibers 
are  considered  to  have  a  slender  and  uniform  cross 
section  that  only  plane  waves  propagate  uniaxially. 
Furthermore,  crossovers  are  modeled  as  hinged  connec¬ 
tions;  then,  slippage  is  assumed  to  be  negligible. 

Mathematical  Formulation 


Considering  a  typical  crossover  in  the  panel/ 
as  shown  in  Figure  9,  the  impulse-momentum  balance 


k^  -Az(j+l,k+li  t) 


Figure  9*  Free-body  diagram  of  forces  acting  at  a  fabric  crossover 
point,  showing  the  influence  of  the  four  fiber  elements 
meeting  there  and  the  elastic  resistive  force  provided  by 
the  fabric  backing. 


during  time  dt  may  be  written  as 
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'Wl  • 


Akt 


(19) 


where  dv  is  an  incremental  velocity,  T  is  the  resultant 
tensile  force,  and  Am  is  the  lumped  mass  of  a 
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fabric  element.  This  relation  provides  a  means  of 
calculating  current  velocity  from  field  variables  in 
previous  time  increment.  For  instance,  at  node  (j+1, 
k+1) ,  the  velocity  at  time  tm+^  may  be  expressed  in  a 
finite  difference  form  as 


•yvt+A  -yyt  ,  A4- 

cr  .  =  AT  ,  +.T  .  ~ L_ 


(20) 


where  £>  is  the  linear  fiber  density,  ,  the 

length  of  an  orthogonal  fiber  element,  £  ,  a  numer¬ 
ical  factor  associated  with  the  crimping  and  wave  of 
the  fabric  structure,  and  the  tensile  force  T  is  given 

A-» 
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(21) 


in  which  is  a  backup  spring  constant  and  £  is  the 

displacement  vector.  T^g  and  Tg^  are  tensile  vectors 

in  the  deformed  orthogonal  fibers  running  through  the 

crossover  as  shown  in  Figure  9.  The  Lagrangian 
*  The  ^underline  is  used  to  denote  a  vector  quantity. 
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coordinates  of  the  node  are  then  evaluated  by 


'vKa-A 

~  k+A 


r 

/N* 


-VYl 

3+1,  k.*A 


•wu-tA 

+AT. 


U*A 


(22) 


The  up-to-date  strain  defined  as  e(j+l,  k+1?  m+1)  = 
^e^g,  eQ^  may  now  be  determined  from  a  continuity 
condition; 
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The  total  kinetic  energy  loss  of  the  missile 


during  penetration  can  be  computed  from 


"v-iV'f-'O 


(28) 


where  V  is  the  residual  velocity  of  the  missile  upon 
penetration.  The  fabric  energy  absorption  and  parti¬ 
tion  are  of  great  interest,  since  they  provide  a  means 
of  evaluating  the  performance  of  the  material.  The 
kinetic  and  strain  energy,  AE  (t)  and  AEf(t), 
obtained  by  the  panel  during  impact  may  be  calculated 
by 


AE ;  cv>  ■=  \  \  4:*  a-w.  I  v-  *  \  JAcl-t  (29) 
At.’ct'i.j  \lT-f  Ca/zeA  cJActt  (30> 

»  A 

Where  A  is  the  area  of  the  fabric.  The  above  formula¬ 
tions  have  been  coded  in  FORTRAN,  and  the  computation 
algorithm  proceeds  from  one  node  to  the  next  along  a 
wave  front  propagating  through  the  fabric.  Due  to 


geometric  symmetry,  only  half  of  one  quadrant  is  con- 


sidered.  At  each  time  increment,  the  code  begins  at 
the  origin  (the  impact  point)  and  works  outward  until 
the  wavefront  is  reached.  The  space  progression  con¬ 
sists  of  a  series  of  passes  along  lines  diagonal  to 
the  orthogonal  fibers,  indicated  by  dashes  in  Figure 
10.  The  algorithm  begins  at  the  x-axis  (a  subroutine 
is  employed  to  handle  the  slightly  different  condi¬ 
tions  along  this  symmetry  boundary) ,  and  progresses 
along  the  diagonal  until  reaching  the  end  of  the 
octant.  For  reasons  of  stability,  the  rate  at  which 
the  code  progresses  outward  from  the  impact  point  is 
related  to  the  fabric  wavespeed?  these  stability  con¬ 
siderations  will  be  discussed  in  a  later  section. 
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Figure  10.  Propagation  scheme  for  the  iterative  wave  propagation 
algorithm.  og 


• '  •  fflSSffiSaB 


jsrpaggyxv^?* 


The  value  of  c?r  is  not  known  prior  to  the  analysis ,  but 
a  preliminary  study  by  Roylance  [15]  indicates  that  it 
is  a  fixed  fraction  of  the  wave  velocity  ir  a  single 
fiber,  cf,  i.e. 

Cp  =  /  OL  (34) 

where  CL  is  a  numerical  factor.  It  generally  has  ■ 
value  greater  than  unity,  which  may  be  attributed  «_o 
the  effective  increase  on  lineal  density  caused  by 
fiber  crossovers.  In  a  square-woven  fabric,  the  lineal 
density  of  a  fiber  along  which  a  wave  is  propagating 
is  effectively  doubled.  This  retards  the  wave  velocity 
according  to  Equation  2  by  a  factor  of  CL  =  "yJlT.  The 
stability  condition  is  then  of  the  form 

At  /tdL&S.')  <1  (35) 

In  the  current  analysis,  Aft  is  obtained  by  physical 
considerations;  therefore  At  is  constrained  by  Equation 
(35).  The  parameter  OL  has  its  optimum  value  obtained 
from  stable  solutions  and  is  defined  as 

DC  =  C  /A-t  (36) 
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In  the  absence  of  a  complete  theory  for  finding 
the  exact  value  of  this  parameter,  it  is  a  common 
practice  in  the  use  of  numerical  methods  to  have  the 
computer  determine  it.  The  optimum  5  may  be  obtained 
by  changing  its  value  continuously  in  test  cases  until 
a  known  solution  is  matched  and  a  variation  of  the 

quantity  will  not  yield  any  appreciable  difference  in 
the  results.  Unfortunately,  there  is  no  existing  solu¬ 
tion  for  this  problem  and  one  must  resort  to  numerical 
tests  of  smoothness  or  of  conservation  laws  for  this 
purpose.  A  convenient  measurement  of  stability  and 
convergence  is  to  study  the  rate  of  energy  conserva¬ 
tion  of  the  system.  In  this  study,  an  energy  discrep¬ 
ancy  parameter  is  introduced  for  the  purpose,  and  is 
defined  as 

^  /E>  ^  1  (37) 

where  Ef  is  the  total  energy  absorption  by  the  fabric, 

and  E  is  the  projectile  energy  loss  defined  previously. 

P 

Figures  11  and  12  give  illustrations  of  the  dependence 
of  solution  stability  and  convergence  on  the  values  of 
.*  parameter.  An  optimum  value  *  =V2  is  obtained 
from  these  results,  in  agreement  with  physical  reasoning. 
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Figure  12.  Illustration  that  the  numerical  scheme  converges  to 

accurate  values  with  time,  as  indicated  by  the  energy 
discrepancy  ratio.  Note  that  nonoptinrum  values  of 
the  stability  ratio  (n  in  this  figure)  lead  to  diver 
gence  at  longer  times. 


Assessment  of  accuracy  of  the  numerical  analysis 


is  somewhat  problematical,  as  no  closed-form  mathematic 


al  analyses  are  available  against  which  to  check  the 


code  results.  Certain  experimental  observations  are 


available,  however,  one  of  which  is  shown  in  Figure  13. 


This  figure  is  a  plot  of  residual  projectile  velocity 


MISSILE  EXIT  VELOCITY,  M/SEC 


Figure  13.  Illustration  that  the  numerical  scheme  predicts  values  of 
final  projectile  velocity  after  penetration  in  agreement 
with  experimental  observation. 


after  penetration  of  a  Kevlar  panel,  as  a  function  of 
initial  velocity.  The  good  agreement  of  the  predicted 
and  observed  results  is  important,  since  it  provides 
some  assurance  that  both  the  transient  response  and 
the  final  fracture  processes  are  being  modeled  reason¬ 
ably.  It  might  also  be  mentioned  that  this  particular 
plot  is  one  which  plays  an  important  role  in  the  design 
process,  so  that  the  ability  to  generate  it  numerically 
without  prior  ballistic  data  or  any  idealizing  assump¬ 
tions  is  of  considerable  practical  importance. 


Another  result  of  the  numerical  calculations  which 
may  be  checked  against  experiment  has  to  do  with  the 
shape  of  the  transverse  deformation  cone,  since  the 
cone  may  be  followed  by  high-speed  photography  during 
the  impact  event.  It  should  be  mentioned  that  this 
photgraphic  evidence  provided  the  initial  impetus  to 
the  development  of  the  present  pin- jointed  fiber 
model,  as  opposed  to  various  membrane  approaches  which 
have  been  attempted  in  the  past.  The  photographs 
clearly  show  a  pyramidal  deformation  cone  which  re¬ 
flects  the  orthogonal  nature  of  the  woven  structure,  as 
opposed  to  the  circular  cones  which  would  be  predicted 
by  axy-symmetric  membrane  analyses.  The  present 
numerical  treatment  predicts  this  pyramidal  shape 

correctly.  A  convenient  indicator  of  deformation  is 
the  size  of  the  cone  at  the  time  of  projectile 
penetration,  as  this  parame  jx  also  reflects  both  tran¬ 
sient  and  fracture  properties  of  the  panel.  Figure  14 
shows  the  predicted  and  observed  cone  size  at  penetra¬ 
tion  for  a  Kevlar  panel,  and  again  it  is  seen  that 
agreement  is  satisfactory. 
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Diameter  at  Penetration,  cm 


Parametric  Materials  Study 


The  utility  of  the  numerical  model  will  be 
illustrated  by  means  of  a  number  of  computer  experiments 
in  which  the  influence  of  fabric  materials  properties 
on  ballistic  resistance  is  assessed.  These  results  help 
validate  the  reliability  of  the  model  in  that  it  can 
be  shown  to  generate  data  in  agreement  with  experimental 
observations.  It  also  provides  a  means  of  illustrating 
certain  phenomena,  such  as  transient  wave  propagation, 
which  are  not  generally  observable  experimentally;  in 
this  regard  one's  intuitive  understanding  of  the  impact 
even  is  improved  considerably. 

Numerical  results  have  been  obtained  for  a 
series  of  four  simulated  orthogonally-woven  square 
panels  203  mm  on  a  side,  impacted  at  zero  obliquity  by 
a  0.22- caliber  projectile  weighing  1.10  gram.  Such  a 

projectile  is  commonly  used  in  experimental  work  to 
simulate  the  effect  of  fragment  impact.  The  edges  of 
the  panels  were  assumed  to' be  clamped,  although  penetra¬ 
tion  generally  occurred  before  the  arrival  of  stress 
waves  at  the  clamps;  the  nature  of  the  edge  boundary 
conditions  is  therefore  relatively  unimportant.  Rather 
than  perform  straightforward  parametric  tests  in  which 
one  variable,  such  as  fiber  modulus,  is  varied  while 
others  are  held  constant,  it  was  decided  to  simulate 


a  series  of  actual  fabrics  for  which  input  data  was 
available  either  from  the  weaver  or  from  laboratory 
measurement.  The  computer  results  can  thus  be  compared 
directly  with  laboratory  ballistic  tests,  although  in 
general,  more  than  one  variable  is  changed  in  each 
simulation.  In  particular,  the  fabric  panel  weight 
varies  slightly  for  each  material  type,  although  this 
effect  was  expected  to  be  small  relative  to  the  large 
change  resulting  from  the  markedly  different  fiber 
moduli . 

For  the  purpose  of  these  parametric  tests,  only 
very  simple  constitutive  and  fracture  models  were 
employed.  Although  more  realistic  models  are  available 
as  described  elsewhere  in  this  report,  the  numerical 
data  necessary  for  input  into  these  models  are  generally 

not  available.  For  this  reason  the  fiber  stiffness  was 
set  to  a  constant  value  obtained  from  handbook  quasi¬ 
static  stress-strain  data,  and  the  failure  criterion 
was  a  simple  maximum-breaking-strain  check,  where  the 
maximum  allowed  was  also  taken  from  quasistatic  tensile 
test  results.  In  spite  of  these  questionable  choices, 
the  results  of  the  computer  simulations  are  of  con¬ 
siderable  interest. 

The  data  for  the  four  fabric  types  are  shown 

below: 
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GEOMETRIC  AND  MATERIALS  PROPERTIES  USED  IN 


FABRIC  STUDIES 


Tensile 

modulus,gpd 

80 

550 

990 

2650 

Fracture 
Strain,  % 

14,0 

4.0 

2.2 

1.1 

Fabric 

Mass,  gm 

19,53 

17,38 

25.75 

27.09 

Yarn  denier[05Q 

1167 

1485 

1500 

Ycm/cm 

17 

16 

16 

16 

Figure 

15  shows 

typical 

computer  predictions  of 

strain  wave 

profiles  obtained 

at  various 

times  after 

400  m/sec  impact  on  the  various  fabrics.  Unlike  impact 
on  a  single  fiber,  in  which  a  constant  level  of  strain 
is  propagated  outward  from  the  impact  point,  the 
array  of  fiber  crossover  junctions  around  the  impact 
point  in  a  fabric  serves  to  reflect  a  portion  of  the 
outward-propagating  wave  back  toward  the  impact  point. 
As  a  result,  the  strain  is  always  greatest  at  the  point 
of  impact,  and  grows  continuously  with  time  (unless 
the  projectile  is  slowed  appreciably  by  the  panel) . 

Both  the  level  of  strain  and  the  rate  of  propagation 


Nylon 


Kevlar  29  Kevlar  49  Graphite 


are  governed  by  the  fiber  modulus  and  density.  As 
shown  in  the  figure,  graphite  fibers  have  the  highest 
modulus  of  the  four  materials,  and  thus  propagate  the 
lowest  level  of  strain  at  the  highest  rate.  As  the 
modulus  is  decreased,  the  strain  level  is  increased  and 
the  wavespeed  is  decreased. 

j  Single  Layer,  8"  x  8" 

Vp  =  400  m/sec 


0  Nylon,  11.59  microseconds 


Kevlar  29,  5. 80  microseconds 

0  9. 02  microseconds 

\  0  12.  25  microseconds 


Graphite,  5.87 


Distance  from  Impact,  cm 

Figure  15.  Distribution  of  strain  along  orthogonal  fibers  passing 
through  the  impact  point.  Curves  are  drawn  for  various 
fabric  types,  at  various  times  after  a  400  m/sec  impact 
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The  strain  history  in  the  fabric  is  directly 
related  to  missile  striking  velocity  as  indicated  in 
Figure  16.  When  the  velocity  is  greater  than  a 
critical  impact  velocity  Vcr,  strain  at  the  impact  point 
continuously  rises  until  penetration,  duo  to  the-  con¬ 
tinual  arrival  of  wavelets  reflected  from  crossovers 
and  boundaries.  In  contrast,  if  the  velocity  is 
smaller  than  VQrr  the  impact  strain  develops  to  a 
'  level  below  the  breaking  strain  and  remains  relatively 

constant  for  the  rest  of  the  dynamic  process.  Here 
the  effect  of  unloading  due  to  projectile  deceleration 
is  able  to  balance  the  increase  of  strain  due  to  wave 
reflection.  The  wavy  strain  history  in  the  figure  may 
be  caused  by  the  dispersion  and  interaction  c-f  the 
traveling  and  reflected  wavelets,  and  perhaps,  some 
numerical  fluctuation. 


TIME  AFTER  IMPACT,  fi t«C 

Figure  16.  Effect  of  initial  projectile  velocity  on  the  development 
of  strain  at  the  point  of  impact  for  nylon  fabric. 


FRACTIONAL  PROJECTILE  VELOCITY 


Penetration  dynamics  can  also  be  illustrated  by 
the  missile  deceleration  as  shown  in  Figure  17,  where 
reductions  o*  missile  velocity  by  various  fabric 
materials  are  given.  Note  that  the  ability  of  the 
various  fabrics  to  decelerate  the  projectile  increases 
monotonically  with  the  fiber  modulus. 


Time  After  Impact,  microseconds 


Figure  IS.  Energy  absorbed  by  a  Kevlar  29  panel  after  a  400  m/sec 
impact,  illustrating  the  partition  of  impact  energy 
into  kinetic  and  strain  energy  in  the  panel. 


Time  After  Impact  (microseconds) 


Figure  19.  Illustration  of  the  relative  ability  of  the  four  fabric 
types  of  absorb  impact  energy.  The  curves  are  termin¬ 
ated  at  the  right  by  projectile  penetration,  as  indicated 
by  a  maximum-breaking  strain  failure  criterion. 


convenient  indicator  of  panel  ballistic  performance, 
and  Figure  19  illustrates  the  relative  energy  absorp¬ 
tion  capabilities  of  the  four  panel  materials  studied. 

It  is  seen  that  the  high  fiber  modulus  of  the 
graphite  panel  leads  to  a  rapid  rate  of  energy  absorp¬ 
tion,  but  that  fracture  occurs  before  the  panel  has 
been  able  to  extract  as  much  of  the  projectile's 
impact  energy  as  the  lower-modulus  fabrics.  Conversely, 
nylon  requires  a  long  time  to  penetration,  but  the 
energy  absorption  rate  is  too  slow  to  lead  to  a  large 
total  energy  absorption.  The  Kevlar  29  panel  exhibits 
the  best  combination  of  energy  absorption  rate  and 
long  time  to  penetration,  and  is  thus  predicted  to  be 
the  superior  ballistic  material  of  the  four  types 
studied. 
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It  should  be  mentioned  that  the  accuracy  of  these 
results  is  limited  by  the  questionable  assumptions 
which  had  to  be  made  due  to  the  present  lack  of  know¬ 
ledge  as  to  dynamic  fiber  properties  which  could  be  used 
as  input  data  for  the  code.  For  these  rather  stiff 
fibers,  the  use  of  a  linear  elastic  constitutive  law 
is  probably  not  a  serious  error;  however,  the  use  of 
maximum-breaking-strain  failure  critrion  is  almost 
certainly  tc  blame  for  some  inconsistencies  in  the 
results.  In  particular,  Kevlar  4  9  is  known  to  be 
essentially  as  good  as,  if  not  superior  to  Kevlar  29 
as  a  ballistic  material.  The  authors  feel  the  shape 
of  the  energy  absorption  curves  in  Figure  19  is 
accurate,  but  that  the  location  of  the  failure  point 
is  poor  in  the  case  of  Kevlar  49.  This  points  out  the 
need  for  more  complete  dynamic  fracture  data  on  these 
fibers,  so  that  more  realistic  models  such  as  that 
described  by  Equation  17  may  be  employed. 

It  is  natural  to  seek  some  simple  relationship 
between  fiber  material  properties  and  fabric  ballistic 
resistance.  The  preceeding  results  lead  one  to  expect 
that  the  most  important  parameter  governing  the  stress 
history  in  the  fabric  before  fracture  is  the  fiber 
modulus.  The  modulus  controls  wavespeed  through  the 
relation  c=  VkE“,  and  thus  the  distance  the  impact  dis¬ 
turbance  will  have  traveled  in  a  given  time.  The 


56 


modulus  also  controls  the  level  of  strain  which  will  be 
generated  by  impact  at  a  given  velocity.  The  relation 
is  not  known  explicitly  for  fabrics,  but  can  be  deter¬ 
mined  by  performing  computer  experiments  using  the 
numerical  code. 

Figure  20  depicts  the  computed  strain  history  at 
the  point  of  impact  for  400  m/sec  impacts  upon  the  four 
model  fabrics.  It  is  clearly  seen  that  an  increase  in 
fiber  modulus  decreases  the  strain  for  a  given  time,  in 
correlation  with  the  same  result  for  single  fibers. 

The  fabric  impact  is  considerably  more  complex  than 
single-fiber  impact,  however;  the  point  of  impact  feels 
not  only  the  continuing  influence  of  the  projectile, 
but  is  also  continually  bombarded  by  wavelets  reflected 
and  diverted  from  adjacent  fiber  crossovers.  The 
situation  is  too  complex  to  permit  simple  generaliza¬ 
tions,  but  the  nonlinear  form  of  the  strain  histories 
for  the  various  fabrics  can  be  taken  to  reflect  the 
influence  of  wave  interactions  occurring  in  a  region 
whose  size  increases  quadratically  with  £irae .  Note 
also  that  the  shape  of  the  strain  histories  varies 
consistently  with  fiber  modulus:  the  time  for  arrival 
of  the  first  peak,  for  instance,  decreases  monotonically 
with  modulus. 
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The  shift  of  the  curves  along  the  abscissa,  however, 
is  less  clear.  The  rate  at  which  the  strain  increases 
at  the  impact  point  is  governed  by  the  complex  inter¬ 
actions  of  waves  traveling  about  within  the  constantly 
expanding  region  of  influence,  and  is  beyond  simple 
visualization.  On  average,  the  time  necessary  for  a 
wave  to  reflect  and  return  eventually  to  the  impact 
point  should  decrease  inversely  with  the  wavespeed, 
i.e.  inversely  with  the  root  of  the  modulus.  However, 
the  size  of  the  region  in  which  stress  waves  are 
traveling  at  any  given  time  also  depends  on  the  wave- 
speed,  and  one  would  expect  that  a  larger  region  of 
influence  would  decrease  the  rate  at  which  reflected 
and  diverted  wavelets  are  able  to  return  to  the  impact 
point. 

It  is  found  that  the  time  after  impact  at  which 
the  first  peak  in  the  strain  occurs  varies  linearly, 
with  good  correlation,  with  the  fourth  root  of  the 
fiber  modulus  (or  the  square  root  of  the  wavespeed) . 
Using  this  observation,  which  is  likely  related  to  the 

geometry  of  the  region  of  influence,  one  can  compensate 

.  0  9  R 

the  abscissal  values  of  Figure  20  by  the  factor  E  . 

The  result  of  the  ordinal  and  abScissal  normalization 

is  shown  in  Figure  21,  where  a  curve  valid  for  all  four 
fabrics  is  developed. 
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Figure  21.  "Master"  curve  for  impact-induced  strain  at  the  point  of 
impact.  Ordinal  values  represent  strain  normalized  on 
the  basis  of  the  strain  which  would  be  generated  in  a 
single  fiber  by  impact  at  the ' same  velocity,  while  abscissal 
values  are  adjusted  by  a  factor  equal  to  the  fourth  root  of 
the  fiber  modulus. 


This  master  curve  represents  an  improved  means 
of  performing  preliminary  armor  design.  Since  the 
normalizing  factors  are  known  once  the  dynamic  modulus 
of  the  fiber  is  specified,  one  can  generate  a  strain 
vs.  time  curve  from  Figure  21  applicable  for  a  particu¬ 
lar  fabric  and  impact  velocity.  The  time  for  rupture 
is  the  time  at  which  the  impact-induced  strain  exceeds 


the  fiber's  dynamic  breaking  strain.  As  in  the  fiber 
case,  we  then  see  that  ballistic  resistance  is  a 
balance  between  high  fiber  modulus  leading  to  high 
wavespeeds  and  lower  strains,  and  fiber  breaking  strain. 

This  approach  is  approximate  in  severax  respects, 
however,  and  is  thus  limited  to  preliminary  design. 

First,  it  is  seen  in  Figure  21  that  perfect  correlation 
among  all  four  test  fabrics  is  not  attained,  the  nylon 
showing  a  deviation  at  high  strain,  similar  deviations 
in  other  fabrics  might  be  observed  as  well.  Second, 
the  curve  of  Figure  21  was  generated  from  computer 
experiments  at  relatively  high  velocity,  so  that  pro¬ 
jectile  slowdown  was  not  an  appreciable  factor.  At 
low  impact  velocities,  the  fabric  is  able  to  decelerate 
the  projectile  and  even  bring  it  to  rest.  The  effect  of 
projectile  slowdown  is  to  generate  unloading  waves  in 
the  fabric  which  travel  simultaneously  with  those 
previously  described.  This  unloading  would  have  a 
strong  influence  on  the  curves  such  as  that  in  Figure 
21,  causing  the  curve  to  pass  through  a  maximum  and 
decrease  thereafter  in  those  cases  in  which  the  fabric 
is  able  to  defeat  the  projectile.  For  these  cases, 
complete  treatment  using  the  numerical  code  would  be 
necessary. 
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III. 


EFFECT  OF  VISCOELASTIC  MATERIALS  RESPONSE 


Viscoelastic  Constitutive  Relations 

In  the  course  of  the  iterative  calculations 
described  earlier,  a  constitutive  material  law  must 
be  evoked  at  each  element  in  order  to  compute  the 
element  tension  from  its  strain  (or  strain  history) . 

One  would  expect  that  a  model  incorporating  viscoelastic 
effects  would  be  necessary  for  proper  simulation  of 
polymeric  materials,  and  in  fact  there  is  considerable 
direct  evidence  [16]  that  relaxation  does  indeed  occur 
in  the  ballistic  time  frame.  This  is  also  to  be 
expected  in  light  of  the  dynamic  mechanical  spectrum  of 
nylon,  for  instance,  in  which  a  beta  relaxation  is 
observed  having  an  apparent  activation  energy  of 

approximately  60  kJ/mole  [17] ;  this  relaxation  is 
calculated  to  occur  in  approximately  five  microseconds 
at  room  temperature. 
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Figure  22.  Wie chert  spring-dashpot  model  for  linear  visco-elastic 
fiber  response. 

A  general  viscoelastic  model  well  suited  for 
computing  tensions  from  pres  :ribed  strains  is  the 
Wiechert  model,  depicted  schematically  in  Figure  22. 
This  model  takes  the  polymer  response  to  be  analogous 
to  that  of  an  array  of  Newtonian  dashpots  and  Hookean 
springs.  The  differential  tension-strain  law  for  the 
jth  arm  of  the  model  is 


where  the  dots  indicate  time  differentiation,  or  is  the 
tensile  stress  and  £  is  the  strain.  Casting  this 
equation  in  finite  difference  form  relative  to  a 
discrete  time  increment  At  and  solving: 


g  )  -cr 


(39) 


where  the  superscripts  t  and  t-1  indicate  values  at  the 
current  and  previous  times  , respectively  ,and  r£  ^  - 

is  a  relaxation  time  for  the  jth  arm.  The  total 
tension  at  time  t  is  the  sum  of  all  the  cr^  plus  the 
tension  in  the  equilibrium  spring  k  : 


(5 


V'  W  - 


1+ 


(40) 


This  tension-strain  calculation  is  performed  at  each 
element  node.  In  addition  to  storing  all  the  k^  and 

j ,  the  computer  must  also  store  the  previous  strain 
and  tension  values  at  each  node. 


The  choice  of  the  kj  and 


j  should  be  such  as 


to  model  the  polymer  viscoelastic  response  in  a  time 
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scale  comparable  to  the  ballistic  event,  which  takes 
place  on  a  microsecond  time  scale.  It  is  of  course 
difficult  to  conduct  such  conventional  tests  as  creep 
or  stress  relaxation  on  this  time  scale,  but  guidance 
as  to  proper  model  parameter  selection  can  be  obtained 
from  dynamic  mechanical  spectra,  using  the  activation 
energies  of  the  appropriate  low-temperature  relaxations 

to  effect  a  temperature-rate  conversion.  For  nylon 
fibers,  for  instance,  one  would  fit  the  Wiechert  model 
to  the  beta  relaxation,  ignoring  the  alpha  and  gamma 
relaxations  as  not  being  appropriate  to  the  ballistic 
time  scale  at  room  temperature. 

Results  for  Single  Fibers 

As  a  means  of  developing  a  proper  context  for 
the  study  of  viscoelastic  response  of  a  woven  textile 
panel,  some  results  obtained  in  an  earlier  study  [18] 
which  used  a  direct  numerical  simulation  of  visco¬ 
elastic  relaxation  in  a  transversely  impacted  single 
fiber  will  be  reviewed  briefly.  The  numerical  approach 
for  this  study  was  identical  to  that  described  for 
fabric  structures,  except  that  it  considered  a  single 
fiber  discretized  as  a  series-connected  assemblage  of 
pin- jointed  finite  elements.  As  in  the  fabric  case, 
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this  treatment  produced  numerical  values  for  the 
position,  veloci  y,  strain,  and  tension  of  each  finite 
element  of  fiber  as  a  function  of  time  after  impact. 

A  variation  of  this  treatment  will  be  described  in  some 
detail  in  Chapter  IV  of  this  report. 


Figures  23  and  24  show  the  distribution  of 
nondimensionalized  strain  and  tension  along  the  fiber 
at  various  times  after  impact,  plotted  against  the 
Lagrangian  fiber  coordinate.  These  distributions  were 
obtained  from  the  Wiechert  model  using  only  a  single 
spring-dashpot  arm  in  parallel  with  the  equilibrium 
spring;  this  *-hree-element  model  is  commonly  known  as 
the  "standard  linear  solid",  or  the  "Zener  solid"  The 
distributions  for  these  two  figures  are  for  a  choice 
of  model  parameters  ke  =  80  gm/den,  k-^  =  20  gm/denf> 
and  rC1  =  50  ^tfsec.  The  values  of  the  ordinates 
have  been  normalized  by  the  strain  or  tension  which 
the  rate-independenc  Smith  theory  predicts  for  a  linear 
elastic  material  at  the  same  impact  velocity. 


Figure  24.  Normalized  tension  distribution  along  fiber 


Tfc.tv-t'i  -To  «*r  L-\h/zJ 


(41) 


£  Cc\=^  «=  e0  ev-p  C-  (42) 

where  +  kg)  is  the  relative  strength  of 

the  viscoelastic  relaxation.  The  wavefront  tension 
magnitude  predicted  by  Equation  41  is  shown  in  Figure 
25,  which  also  serves  to  illustrate  the  numerical 
accuracy  of  the  direct  analysis.  (Here  a  0.15  m 
fiber  was  divided  into  200  finite  elements) .  Some 
numerical  overshoot  is  evident  at  the  discontinuous 
wavefront,  but  the  distribution  extrapolates  to  the 
analytically-predicted  value. 

By  means  of  Laplace  transforms.  Smith  [20]  also 
obtained  approximate  expressions  for  the  strain  and 
tension  distributions  in  a  longitudinally  impacted 
fiber.  These  expressions  predict  that  the  tension  and 
strain  at  the  point  of  impact  will  approach  the 
limiting  values 

TCo^  “T0  =  c?«  “To  (43) 
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Where  the  numerical  coefficients  are  for  Vs  =  0.2. 

At  x  =  0,  the  distributions  in  Figures  25  and  23 
approach  limiting  values  greater  than  Equation  43  for 
tension  and  greater  than  Equation  44  for  strain.  Thus 
stress  relaxation  is  slightly  less  and  creep  slightly 
greater  for  transverse  impact  than  for  longitudinal 
impact;  Smith  [19]  reached  this  same  conclusion  in  his 
work  on  transverse  impact. 


X  (cm) 


Figure  25.  Numerical  values  for  tension  distribution  for  t  -  41.08 
microsec  after  intact. 


Results  for  Woven  Panels 


The  three-element  Wiechert  model  (the  standard 
linear  solid)  used  in  the  previous  section  was  also 
employed  to  examine  the  influence  of  viscoelastic 
relaxation  during  ballistic  impact  of  woven  panels. 

The  model  parameters  were  chosen  to  simulate  ballistic 
nylon:  k^  =  20  gm/den,  kg  =  80  gm/den,  *C  ^  -  5  sec. 
Results  have  been  obtained  for  a  simulated  0.2  m  x  0.2-m 
panel  weighing  19.5  gm,  impacted  with  a  0.22-caliber 
fragment  simulating  projectile  weighing  1.10  gms  at 
various  impact  velocities. 


°0  I  2  3  4  5  6  7° 

DISTANCE  FROM  IMPACT  POINT,  cm 

Figure  26.  Stress  distributions  along  orthogonal  fibers  running  through 
impact  point  for  linear  elastic  and  viscoelastic  fabrics  (t  •• 
30.4  microsec) . 
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Figure  26  shows  the  distribution  of  stress  along 
the  orthogonal  fibers  running  through  the  point  of 
impact  at  t  =  30.4  ^tsec  after  an  impact  at  Vp  =  300 
m/sec.  The  results  for  the  viscoelastic  fabric  are 
compared  with  those  of  an  ideally  elastic  fabric 
having  a  stiffness  equal  to  that  of  the  unrelaxed 
viscoelastic  material  (100  gm/den) .  The  nonuniform 
distributions  along  the  fiber  are  due  as  described 
earlier  to  the  continual  reflection  of  wave  components 
from  fiber  crossovers,  resulting  in  a  maximum  in  stress 
at  the  impact  point.  An  appreciable  difference  in 
stress  levels  between  the  linear  elastic  and  Visco¬ 
elastic  cases  is  observed,  especially  near  the  wave- 
front.  Relaxation  of  the  stress  due  to  the  rate- 
dependent  material  behavior  may  be  expressed  by  the 
ratio  of  the  viscoelastic  and  the  elastic  and  the 
elastic  stress,  T/T  ,  as  shown  in  the  figure.  It  is 
found  that  a  large  amount  of  relaxation  occurs  near 
the  wave  front  while  an  equlibrium  state  of  relaxation 
is  reached  in  the  region  away  from  the  disturbance  front. 
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DISTANCE  FROM  IMPACT  POINT,  cm 

Figure  27.  Distribution  of  strain  along  orthogonal  fibers  running 
through  impact  point. 

The  wave  attenuation  during  this  300  m/sec  impact 
is  also  demonstrated  in  figure  27,  where  relative  strain 
distributions  fc  /&^  are  given  for  various  times 
after  impact.  As  illustrated  in  this  figure,  the 
magnitude  of  the  wavefront  attenuates  significantly  as 
it  propagates  along  the  orthogonal  fibers,  and  the 
strain  at  a  given  position  increases  with  time  after 

impact  from  its  original  values. 
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RELAXATION  AT  IMPACT 


TIME  AFTER  IMPACT,  ft  sec 

Figure  28.  Stress  histories  at  impact  point  for  linear  elastic  and 
viscoelastic  materials. 


The  stress  and  strain  histories  at  the  point  of 
impact  shown  in  Figure  28  give  another  indication  of 
viscoelastic  dissipation  in  the  response  of  the  panel. 
They  increase  continuously  with  time  due  to  the  re¬ 
flection  of  wavelets  from  crossovers;  however,  stress 
relaxation  and  strain  creep  of  the  viscoelastic  material 
occur  simultaneously  with  this  general  increase.  The 
viscoelastic  stress  at  a  given  time  is  smaller  than  the 
elastic  case  as  shown  in  the  figure.  The  relative 
relaxation  at  the  point  of  impact,  denoted  by  the  ratio 
of  the  viscoelastic  and  elastic  case,  develops  gradually 
and  reaches  a  steady  state  at  long  times.  The  relaxa¬ 
tion  histories  for  different  missile  striking  velocities 


RELAXATION  AT  IMPA 


are  given  in  Figure  29.  The  similarity  of  their  magni 
tudes  is  a  manifestation  of  the  linear  material  res¬ 
ponse.  Again,  these  stresses  approach  an  asymptotic 
equilibrium  state  at  times  longer  than  the  characteris 
tic  material  relaxation  time. 


0.8 


Vp  =  400  m*ac  *  Vp=300  m/aec  *  *  Vp  «  250  m/sec 
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TIME  AFTER  IMPACT,  psec 

Figure  29.  Stress  relaxation  at  impact  point  for  various  impact 
velocities. 

Nonlinear  Viscoelastic  Response 

The  use  of  the  Wiechert  model  as  described  in 
the  previous  sections  is  sufficient  to  illustrate  the 
most  important  features  of  rate-dependent  ballistic 
materials.  However,  most  materials  do  not  meet  the 


rigid  requirements  of  linearity  necessary  for  a  truly 
rigorous  application  of  concepts  of  linear  viscoelasti¬ 
city.  For  more  detailed  simulations  of  fabric  and  fiber 
response,  one  must  turn  to  general  constitutive  rela¬ 
tionships  which  more  accurately  model  the  behavior  of 
these  materials.  Unfortunately,  there  does  not  exist  a 
general  concensus  of  the  most  realistic  means  of 
achieving  this  goal.  At  present,  the  subject  of  non¬ 
linear  viscoelasticity  is  being  pursued  actively  by 
several  groups,  several  of  which  are  employing  highly 
divergent  approaches  to  this  problem.  It  is  not 
possible  here  to  select  any  single  approach  as  having 
significantly  greater  merit  than  certain  others. 

However,  the  ease  with  which  various  constitutive 
laws  may  be  incorporated  into  the  direct  analysis  scheme 
makes  it  possible  to  assess  relatively  easily  the  in¬ 
fluence  of  various  assumed  material  models.  As  an 
illustration  of  this  capability,  some  results  using 

Eyring's  model  of  thermally-activated  nonlinear 
viscoelasticity  will  be  presented. 

The  computationally-convenient  Wiechert  model 
can  be  extended  to  include  the  effect  of  material  non¬ 
linearity  by  rendering  the  springs  and/or  the  dashpots 
nonlinear.  If, for  instance  ,one  uses  a  power-law  spring 
and  a  nonlinear  Eyring  dashpot  [21]  ,  defined  as 
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then  the  finite-difference  equation  relating  tensions 
and  strains  in  the  jth  arm  of  the  model  is: 
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A  relation  such  as  this  requires  an  iterative  numerical 
solution  for  CTj  at  each  element  and  at  each  time 
step;  the  computer  effort  is  increased  but  the  princi¬ 
ples  of  the  impact  algorithm  are  straightforward.  The 
principal  obstacle  to  the  use  of  nonlinear  models  in 
the  direct  analysis  is  not  the  incorporation  of  the 
models  into  the  computational  scheme,  but  rather  the 
determination  of  the  material  parameters  (the  b’s , 
k 1  s  ,A '  s  ,and  oi.'s  in  Equation  46)  applicable  to  the 

microsecond  time  scale  of  polymer  relaxations. 

To  illustrate  the  effect  of  nonlinear  constitutive 
models  on  panel  ballistic  response,  a  series  of  computer 
experiments  was  performed  on  three  different  simulations 
of  nylon  fabric:  one  using  only  linear  elastic  response 
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(only  the  equlibrium  spring  in  the  Wiechert  model) ,  one 
using  the  standard  linear  solid  model  for  linear  visco¬ 
elastic  response  (the  equilibrium  spring  plus  one  spring- 
dashpot  arm) ,  and  the  last  being  a  standard  linear 
solid  but  with  the  dashpot  made  a  nonlinear  Eyring 
element.  The  intial  modulus  was  taken  as  100  gm/den, 
the  relaxed  modulus  as  80  gm/den,  and  the  relaxation 
time  for  the  standard  linear  solid  as  five  ^K.sec  (the 
same  as  in  the  previous  section) .  The  concept  of  rela¬ 
xation  time  (the  time  to  complete  63.2%  of  the  total 
relaxation)  has  no  meaning  for  the  nonlinear  element, 
since  the  rate  of  relaxation  changes  nonlinearly  with 
the  stress.  Lacking  any  experimental  data  in  this  time 
scale,  the  A  and  oi  were  arbitrarily  chosen  so  as  to 

cause  relaxation  in  approximately  the  same  time  scale 

as  the  standard  linear  solid.  A  and  CL  were  set  at 
3  —1 

10  sec  and  0.7  den/gm, respectively.  The  nonlinear 
constitutive  equation  was  solved  at  each  element  using 
Muller ’s method  [22],  which  increased  the  computation 
time  relative  to  that  of  the  standard  linear  solid  by 
roughly  one-third. 
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DISTANCE  FROM  IMPACT  POINT,  cm 

Figure  30 •  Stress  distributions  along  orthogonal  fibers  running 
through  impact,  point  for  linear  elastic  and  nonlinear 
viscoelastic  fabrics. 

One  means  of  comparing  the  various  constitutive 
models  is  in  terms  of  the  distributions  of  strain  along 
the  orthogonal  fibers  running  through  the  impact  point, 
at  various  times  after  impact.  Figure  30  shows  the 
distributions  for  the  elastic  and  nonlinear  viscoelastic 
materials  20  ^ttsec  after  a  300-m/sec  impact.  Also 
shown  is  the  ratio  between  the  nonlinear  and  the 
elastic  cases.  This  ratio  is  a  measure  of  the  stress 
relaxation  in  the  fabric;  it  is  greatest  at  the  wave- 
front,  as  the  large  gradient  of  strain  there  produces  a 
similarly  large  rate  of  relaxation. 
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DISTANCE  FROM  IMPACT  POINT,  cm 

Figure  31  •  Comparisor.  of  stress  relaxation  in  linear  and  nonlinear 
viscoelastic  fabrics. 

The  nonlinear  stress  relaxation  is  compared  in 
Figures  31  to  that  produced  by  the  standard  linear 
solid.  Although  the  equilibrium  value  far  from  the 
wavefront  is  approximately  the  same  in  both  cases, 
strong  differences  are  evident  near  the  wavefront.  These 
are  due  to  the  relatively  more  rapid  response  to  higher 

strain  gradients  in  the  nonl'-.ear  material.  Another 
indicator  of  viscoelastic  fabric  response  is  the  stress 
at  the  point  of  impact.  The  stress  and  strain  at  the 
impact  point  increase  with  time-  due  to  the  continual 
arrival  there  of  wavelets  reflected  from  fiber  cross¬ 
overs,  but  in  viscoelastic  materials  both  stress  relax- 
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Figure  32.  Stress  histories  at  impact  point  for  linear  elastic, 

linear  -viscoelastic,  and  nonlinear  viscoelastic  fabrics. 
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IV.  NUMERICAL  ANALYSIS  OF  WAVE  PROPAGATION 


IN  TWO  CROSSED  FIBERS 

Introduction 

This  chapter  describes  study  of  the  dynamics  of  a 
special  but  highly  important  physical  system:  that  of 
two  fibers,  one  having  been  transversely  impacted  at 

zero  obliquity  by  a  high-speed  projectile,  and  the 
other  crossing  the  first  perpendicularly  at  some  . 
distance  from  the  impact  point.  This  system  is  germane 
to  the  understanding  of  impact  and  wave  propagation 
phenomena  in  woven  textile  panels  used  for  ballistic 
protection.  The  wave  propagation  phenomena  occurring  at 
the  fiber  crossover  have  a  strong  influence  on  the 
response  of  a  woven  panel  to  impact,  since  these  panels 
typically  have  on  the  order  of  forty  crossovers  per 
inch.  The  nature  of  these  crossover  interactions  may 
be  one  of  the  factors  causing  what  appears  to  be  an 
excellent  fiber  in  single-fiber  ballistic  tests  to 
exhibit  less  ballistic  protection  when  woven  into  a 
textile  panel  structure  than  a  nominally  inferior  fiber. 
As  mentioned  earlier,  this  situation  obtains  in  the 
case  of  the  Kevlar  ballistic  protectio:.  vests  now 
being  used  by  military  and  police  personnel:  the  Kevlar 
vest  outperforms  the  older  nylon  vest,  in  spite  of 
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nylon's  having  a  higher  transverse  critical  velocity. 


The  inability  to  predict  vest  performance  from  single¬ 


fiber  test  data  is  a  matter  of  considerable  concern  to 


the  armor  design  community,  and  this  study  of  fiber 


crossover  dynamics  was  begun  to  clarify  this  situation 


Method  of  Solution 


*  Crossover 


Figure  33*  Schematic  of  model  for  numerical  analysis  of  two  crossed 
fibers . 
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System  Idealization.  The  system  of  two  crossed 
fibers  is  modeled  as  in  Figure  33,  where  the  origin  of 
coordinates  is  placed  at  the  midpoint  of  the  clamped 
primary  fiber,  which  extends  along  the  x-axis.  The 
projectile  moves  along  the  y-axis  only,  and  impacts  the 
primary  fiber  at  the  origin.  From  symmetry,  only  half 
the  primary  fiber  need  be  considered.  The  secondary 
fiber  extends  along  the  z-axis  and  intersects  the 
primary  fiber  at  some  a- -  .^rary  distance  from  the 
origin.  At  the  crossover  point,  the  secondary  fiber  is 
assumed  to  follow  the  motion  of  the  primary  fiber  in 
the  direction  perpendicular  to  the  primary  fiber  (in 
the  x-y  plane) ,  but  is  allowed  some  measure  of  slip  in 
the  direction  parallel  to  the  primary  fiber.  Motion  of 
the  primary  fiber  is  assumed  to  occur  in  the  x-y  plane 
only,  while  the  secondary  fiber  may  move  in  all  three 
directions. 


Figure  34*  Discrete  element  of  fiber. 


I 
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To  proceed,  the  fibers  are  discretized  as  a  series 
of  n  pin-jointed  finite  elements  of  equal  length  as 
shown  in  Figure  34.  The  masses  of  the  elements  are 
taken  to  be  lumped  at  the  nodal  end  points  of  these 
elements,  and  at  these  nodes  are  defined  vector  co¬ 
ordinates  x,-  ,  velocity  v.  ,  and  tension  T.  .  The 

<*->i 

scalar  strain  £  ^  at  each  element  will  be  computed 
from  the  coordinates  of  the  nodes  at  either  end  of 
the  element.  The  tension  T.  has  the  same  direction 
as  the  element  itself  (approximating  the  element's 
assumed  inability  to  support  a  bending  moment) ,  while 
v-  is  not  constrained  in  direction.  These  elements  are 
now  described  as  in  the  fabric  analysis  by  simple 
governing  equations:  impulse -momentum  balance,  strain- 
displacement  relation,  constitutive  relation,  etc. 

These  relations  are  cast  as  a  recursive  algorithm  for 
proceeding  from  one  element  to  the  next  along  the 
fiber  length,  and  then  repeating  the  process  at  a  new 
increment  of  time.  The  computer  solution  thus  is  re¬ 
ferenced  to  a  Lagrangian  frame  of  reference  attached 
to  and  extending  with  the  fiber,  which  effectively 
reduces  the  problem  to  one  dimension. 

Momentum  Balance.  A  consideration  of  impulse- 

st 

momentum  balance  at  the  i  +  1  node  provides  a  means 
of  computing  the  current  velocity  at  that  node  in  terms 
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Of  its  velocity  in  the  previous  time  lnorement  an<J  thfi 
tensions  acting  on  it  during  that  time  increment.  (I„ 
the  folio  ing,  subscripts  on  a  variable  refer  to  the 
which  it  is  defined , while  superscripts  t  and 


t  1  refer  to  values  at  the  current  and  previous  times 
respectively) .  The  impulse-momentum  balance  can  be 
written  in  finite  difference  form  as 


ST  ~  Ay*. 

At 


(47) 
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Letting  A  &m/  At,  a  fixed  parameter,  equation 
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may  be  solved  for  vfc. 
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This  vector  expression  can  be  written  in  scalar  form 

by  reference  to  the  inclination  angles  P.  and  «. 

°f  the  and  Ti+1  vectors  respectively: 
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Then  the  x  and  y  components  of  velocity  are: 
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where  T  =  1  tJ  is  the  tension  magnitude.  The  boundary 
conditions  are  easily  incorporated  into-  the  impulse- 
momentum  balance:  at  the  first  mode,  the  velocity  is 
set  equal  to  the  current  projectile  velocity  (v^  =  V*") , 

and  at  the  clamp  the  velocity  is  set  to  zero  (vt  =  0) . 

**n 


Strain-displacement  Relation.  Having  computed 
the  velocities  at  the  i^*1  and  iHs^  nodes,  the  strain  in 
the  element  between  these  nodes  is  computed’  as 
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polynomial  and  exponential  strain  hardening) ,  linear 
viscoelasticity  (Wiecherc  model) ,  and  nonlinear  visco¬ 
elasticity  (Eyring  model) . 


Computation  of  New  Projectile  Velocity.  The 
algorithm  described  above  proceeds  from  one  element  to 
the  next  along  the  length  of  fiber,  and  i"  started  by 
imposing  the  initial  projectile  velocity  on  the  first 
node.  At  the  end  of  the  first  time  increment,  a  strain 
will  have  developed  in  the  first  element  due  to  the 
velocity  difference  between  the  first  and  second  nodes. 
(Initially,  all  velocities,  tensions, and  strains  are 
set  to  zero.)  This  strain  produces  a  tension  as  cal¬ 
culated  from  the  constitutive  relation,  and  this  tension 
produces  a  velocity  in  the  second  node  beginning  at 
the  next  time  increment. 

After  each  time  increment,  at  the  completion  of 
the  lengthwise  recursive  calculations,  a  new  projectile 
velocity  v^  can  be  computed  by  means  of  a  momentum 

r 

balance  using  the  tension  at  the  first  node: 
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where  M  is  the  projectile  mass,  T  is  the  component  of 
P  y 

fiber  tension  in  the  projectile  travel  direction,  and 
the  factor  2  accounts  for  the  other  half  of  the  primary 
fiber  extending  in  the  -x  direction.  T  is: 
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Crossover  Fiber  Calculations.  Computation  of 
field  variables  along  the  secondary  fiber  proceeds  in  a 
lengthwise  manner  exactly  as  described  above,  although 
the  vector  resolutions  become  slightly  more  complicated 
due  to  the  motion  in  three  rather  than  two  dimensions. 
The  secondary  algorithm  is  started  by  imposing  on  the 
first  node  of  the  crossover  fiber  the  velocity  imparted 
to  it  by  the  primary  fiber.  As  stated  earlier,  the 
secondary  fiber  is  allowed  a  measure  of  slip  along  the 
primary  fiber  but  is  constrained  to  follow  it  in  the 
direction  normal  to  the  primary  fiber.  Denoting  the 
node  on  the  primary  fiber  nearest  the  crossover  point  as 
ix,  the  velocity  of  this  node  resolved  in  directions 
parallel  and  perpendicular  to  the  primary  fiber  there 
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In  equation  62,  the  notation  of  the  form  yl  or  y2 
indicates  field  variables  for  the  primary  and  secondary 
fibers  Respectively.  The  velocity  imposed  on  the  first 
(crossover)  node  of  the  secondary  fiber  is: 
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where 


Qi  is  a  slide  factor  which  permits  no  sliding 
s 

when  set  to  1  and  unrestrained  sliding  when  set  to 
zero. 

The  crossover  node  ix  will  change  with  time  if  the 

tangential  slip  along  the  primary  fiber  is  sufficient. 

After  each  time  increment,  a  new  position  of  the  first 

node  on  the  secondary  fiber  is  computed,  and  ix  is 

assigned  to  the  nearest  node  on  the  primary  fiber. 

The  momentum-balance  calculation  of  v. ..  m  the 

-*'1+1 

,  ,  ,  cf  . 

primary  fiber  must  be  modified  when  the  l+l  -  node  is 
also  the  crossover  node  ix,  since  the  secondary  fiber 
applies  its  own  tension  to  that  node.  Denote  the  direc¬ 
tion  angles  of  the  primary  and  secondary  fibers  at  that 

node  as  <H, » to*  and  to  v4>z  *  ,respect- 

^  *  a  Tl 

ivesly  (  to,  =  0) .  Then  the  components  of  tension 
applied  by  the  secondary  fiber,  resolved  along  directions 
parallel  and  perpendicular  to  the  primary  fiber,  are: 
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pava.  1  L  %  X  Sr  1 


pava. 

4: 


Tt?1  «"TZL  co* 4>1  -  eos<to. eos^M.  1  (67) 

petp  H  \  %  *  V 
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where  the  direction  cosines  are  computed  from  the 

current  nodal  coordinates.  The  primary  fiber  is  allowed 

to  feel  the  impulse  of  the  perpendicular  component 

fully,  but  the  parallel  component  is  reduced  by  the 

slide  variable  6L  .  The  usual  computation  of  v*" 

s  ~i+l 

is  the  adjusted  as 

4r  4;  r  ^ 

'■'^-1  -  .  'm~'  will ,  +  Z.A\^”Z-  Cos  4^1  <u 

1*1  %ft\  L  pevp  ^ 


■t  t  r  t 

.  AT.  tZA  TZ1 

1.+A  f.+A  L  P 


4-  cL^Tz_  CX«. 


p«*r 


where  the 


+  ^"TZ  ^  cos 

symbol  indicates  a  computer  replacement 


operation;  i.e.  the  additional  impulse  from  the 
secondary  fiber  is  added  to  that  already  computed  from 
equations  63  and  64 . 


Stability,  Accuracy  ,and  Efficiency.  Criteria  for 
stability  and  accuracy  of  the  above  method  are  related 
as  in  the  fabric  case  to  the  theory  of  charactaristics 
for  hyperbolic  systems  of  partial  differential  equations 
and  are  similar  to  those  for  finite-difference  solutions 
of  wave  propagation  problems  [14]  .  Given  a  wave 
equation  of  che  form 
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(70) 


which  is  to  be  solved  by  approximating  ^t  and 
by  finite  differences  At  and  A*»  a  stability 
ratio  can  be  defined  as 


(71) 


The  finite  difference  scheme  is  stable  and  accurate  for 
6i  =1,  stable  but  increasingly  inaccurate  for  au.  1 
and  unstable  for  >  1 .  The  choices  for  Ax  and 

At  are  thus  not  independent,  but  are  related  by  the 
wavespeed  for  the  choice  of  0L  =  1 . 

In  the  direct  analysis  of  the  fibers  described 
above,  this  stability  criterion  is  equivalent  to  adjust 
ing  the  rate  of  march  of  the  computer  solution  along 
the  fiber  to  match  the  rate  of  propagation  of  the 
strain  wave.  Conceptually,  this  requirement  is  related 
to  the  necessity  of  programming  the  finite  governing 
equations  so  as  to  model  the  actual  continuous  dynamic 
process  as  accurately  as  possible.  If  a  major  disturb- 
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ance  -  such  as  the  passage  of  a  strain  wave  with  its 
accompanying  energy  input  -  takes  place  in  a  finite 
element  which  is  not  considered  explicitly  in  the 
computational  scheme,  divergent  numerical  results  are 
very  likely. 

Once  a  stable  computational  scheme  has  been  • 
developed,  one  usually  attempts  to  increase  its  accuracy 
to  whatever  limit  is  desired  by  decreasing  the  size  of 
the  elements;  l.e.  oy  increasing  the  number  of  nodes. 
Since  for  6L  =  1  a  decrease  in  Ax  requires  a  corres¬ 
ponding  decrease  in  t,  the  computation  time  -  and 

therefore  the  expense  -  required  for  analysis  of  a 
given  impact  event  increases  as  the  square  of  the 
number  of  nodes.  The  element  size  is  therefore  chosen 
so  as  to  balance  the  conflicting  requirements  of 
economy  and  accuracy.  As  an  example  of  computation 
time,  the  CPU  requirement  for  the  IBM  370/168  system 
was  0.168  minutes  for  a  problem  in  which  the  strain  wave 
propagated  0.2  m  along  the  primary  fiber  and  0.1  along 
the  secondary  fiber,  with  a  length  increment  of  2.0  mm. 
As  a  means  of  improving  code  efficiency,  the  program 
employs  logical  flags  which  terminate  the  length  loop 
computation  when  the  computer  passes  the  point  along  the 
fiber  length  corresponding  to  the  wavefront. 
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Accuracy  assessment  for  the  case  of  two  crossed 
fibers  is  difficult,  since  no  experimental  or  closed- 
form  mathematical  analysis  of  this  problem  is  available, 
but  some  assurance  of  accuracy  is  derived  from  computer 
runs  in  which  the  secondary  fiber  is  placed  at  the 
origin  (the  impact  point).  In  this  case,  response  of 
both  the  primary  and  secondary  fibers  is  found  to  be 
that  predicted  by  independent  analyses.  Data  such  as 
that  previously  presented  in  Figure  23  is  obtained  along 
the  primary  and  secondary  fibers.  In  certain  cases  to 
be  discussed  below,  the  numerical  overshoot  and 
oscillation  observed  near  wavefronts  cause  problems  in 
interpretation  of  results.  These  oscillations  are  a 
result  of  the  inability  of  the  discrete  difference 
equations  to  model  discontinuities  accurately.  Although 
the  method  is  conditionally  stable  and  the  oscillations 
are  damped  out  away  from  the  discontinuity,  problems 
of  interpretation  remain  near  the  discontinuity.  The 
oscillation  at  the  wavefront  is  diminished,  by  the 
material  viscosity,  and  in  some  cases  an  "artificial" 
viscosity  may  be  included  solely  for  the  purpose  of 
smoothing  the  numerical  results. 


The  figure  shows  the  distribution  of  strain  in  each 
fiber  28.7  ^csec  after  impact  at  400  m/sec ,  where  the 
abscissa  measures  the  distance  along  the  secondary 
fiber  from  the  crossover  point.  The  dotted  line  at 
strain  =  1.45%  depicts  the  level  of  strain  which  would 
be  generated  in  a  single  fiber  at  this  impact  velocity. 

In  this  example  no  viscosity  has  been  included,  and  the 
large  overshoot  at  the  wavefront  causes  problems  in 
interpretation  of  results.  In  spite  of  this  oscillatory 
behavior,  however,  an  increase  in  strain  in  the  primary 
fiber  behind  the  crossover  due  to  the  wavelet  reflected 
from  the  crossover  is  evident,  as  is  a  reduction  in  the 
strain  intensity  in  the  region  of  the  primary  fiber 
beyond  the  crossover.  More  easily  measured  is  the  level 
of  strain  intensity  propagated  along  the  secondary  fiber. 

Computer  experiments  were  conducted  on  the 
crossover  system  for  a  range  of  fiber  moduli  and  slide 
factors,  and  graphical  output  similar  to  Figure  35  used 
to  determine  coefficients  of  wave  reflection,  trans¬ 
mittance,  and  diversion.  These  coefficients  are  defined 
as  that  fraction  of  the  outward-propagating  strain  wave 
which  is  reflected  backwards  by  the  crossover,  the 
fraction  which  passes  through  the  crossover  and  continues 
its  outward  propagation,  and  the  fraction  which  is 
diverted  and  begins  propagating  along  the  fiber  passing 
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transversely  through  the  crossover,  hs  a  means  of 
obtaining  these  coefficients  in  spite  of  the  uncertain¬ 
ties  caused  by  the  numerical  fluctuations  near  the 
wavefronts,  the  computer  was  asked  to  determine  the 

away  from  the  oscillation  region.  In  order  to  guarantee 
conservation  of  energy,  the  sum  of  the  sguares  of  the 
above  three  coefficients  should  equal  unity,  this  was 
in  fact  obtained  and  offers  some  assurance  as  to  the 
accuracy  of  the  numerical  values. 
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The  variation  in  the  transmission  and  diversion 
coefficients  with  fiber  modulus  is  shown  in  Figure  36. 
The  coefficient  of  reflection  was  near  1%  over  this 
range  of  moduli,  but  showed  considerable  scatter.  It  is 

seen  that  the  diversion  coefficient  is  of  a  much  larger 
magnitude  than  the  reflection  coefficient,  and  that  it  varies 

more  strongly  with  the  fiber  modulus.  The  major 
portion  of  the  crossover  influence  on  wave  propagation 
is  thus  ascribed  to  diversion  rather  than  reflection. 

This  observation  is  of  significance,  since  an 
approximate  treatment  of  fabric  impact  by  Freeston  and 
Claus  [23]  sought  to  predict  the  increase  of  strain  at 
the  impact  point  by  considering  wave  propagation  along 
a  single  fiber  which  reflects  a  certain  portion  of  the 
outward-propagating  wave  at  a  series  of  discrete  points 
along  its  length.  The  analysis  is  then  reduced  to  a 

bookkeeping  procedure  in  which  one  keeps  track  of 
inward  and  outward-propagating  waves  in  each  of  the 
elements  between  these  reflection  points.  This  scheme 
leads  to  a  very  simple  computer  code  and  one  would  hope 
that  it  could  provide  at  least  approximate  guidance  in 
design. 


Figure  37.  A  comparison  of  the  reflection-only  bounce  model  for  wave 
propagation  in  an  impacted  fabric,  in  comparison  with  the 
fabric  model  of  this  report. 


Unfortunately,  it  appears  that  the  reflection-only 
model  predicts  much  too  high  a  strain  level,  except  at 
times  very  early  in  the  ballistic  event.  Figure  37 
shows  the  impact-point  strain  history  as  predicted  by 
the  reflection-only  model  for  a  400  m/sec  impact  on  a 
Kevlar  29  panel  with  1575  yarns/m  (the  curve  developed 
by  the  "BOUNCE"  code) .  This  prediction  is  compared 
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with  that  of  the  more  general  code  described  earlier 
(the  curve  labeled  "FABRIC").  At  very  short  times,  the 
BOUNCE  code  results  are  likely  superior,  as  they  give 
strains  equal  to  that  developed  in  a  single  transversely- 
impacted  fiber;  the  Fabric  code  shows  a  numerical  lag 
in  the  development  of  strain.  The  two  codes  achieve 
similar  values  at  near  1.5  microseconds,  but  after  this 
the  bounce  code  increases  rapidly  to  unreasonably  high 
values  of  strain  and  thus  predicts  penetration  too 
early.  It  is  interesting,  however,  that  the  value  of 
the  reflection  coefficient  chosen  by  Freeston  and  Claus 
in  order  to  bring  their  model  into  line  with  experiment 
was  very  nearly  that  found  explicitly  in  the  crossover 
study  (0. 01) . 

As  the  slide  factor  dt,  decreases  from  “unity  toward 
zero,  representing  less  fiber-fiber  friction  at  cross¬ 
over  points,  one  would  expect  that  the  reflection  and 
diversion  coefficients  would  approach  zero  and  that 
the  transmission  coefficient  would  approach  unity.  At 
=  0,  there  is  no  coupling  between  the  two  fibers 
(until  the  arrival  of  the  transverse  kink  wave,  which 
generally  occurs  later  than  the  arrival  of  the  longi¬ 
tudinal  wave).  As  seen  in  Figures  38,  39  and  40, 
respectively i for  Kevlar  29  fibers,  this  trend  is 

quantified  by  the  results  of  the  crossover  computations. 
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SLIDE  FACTOR 


Figure  38.  The  influence  of  fiber-fiber  sliding  on  the  fraction 
of  stress  wave  intensity  which  is  reflected  at  fiber 
crossovers,  as  indicated  by  computer  experiments  on' 
Kevlar  29  fibers. 
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Figure  39.  The  influence  of  fiber-fiber  sliding  on  the  extent 
to  which  a  portion  of  the  propagating  stress  wave 
is  diverted  from  the  primary  fiber  to  begin  pro¬ 
pagating  along  the  transverse  secondary  fiber. 


SLIDE  FACTOR 


Figure  2+0 .  The  influence  of  sliding  on  the  extent  of  stress  wave 
intensity  propagated  beyond  fiber  crossovers. 


In  principle,  it  would  be  possible  to  include  the 
effects  of  fiber-fiber  slippage  in  the  two-dimensional 
fabric  code  by  incorporating  the  formulae  of  this 
present  chapter  into  the  general  code.  Such  an  in¬ 
corporation,  however,  would  likely  render  the  fabric 
code  so  much  slower  as  to  be  uneconomic.  In  addition, 
one  has  at  present  no  real  means  of  assessing  the  slide 
parameter  needed  for  the  computation.  It  might  be 
possible,  however,  to  adjust  the  viscosity  of  the 


material  model  artificially  in  order  to  produce  results 
similar  to  those  caused  by  fiber  slippage.  This  method 
would  be  similar  to  that  employed  during  wave  propaga¬ 
tion  calculations  in  geometrically  dispersive  composites 
[24].  If  this  approach  is  pursued  in  the  future,  the 
data  of  Figures  38-40  will  be  useful  in  providing 
guidance  as  to  the  desired  effect.  As  a  final  comment 
on  this  work,  it  may  be  stated  that  the  FABRIC  impact 
code  provides  a  simulation  of  impact  on  textile  fabrics 
which  is  already  of  sufficient  accuracy  that  inclusion 
of  the  fiber  crossover  effects  would  not  be  considered 
necessary,  at  least  in  the  case  of  tightly  woven  panels 
which  do  not  exhibit  extensive  fiber  slippage  during 
impact . 
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CONCLUSIONS 

The  numerical  analyses  described  in  this  report 
offer  a  means  whereby  the  designer  of  personnel  armor 
may  perform  computer-aided  design  and  analysis  of  what 
up  until  now  has  been  an  impossibly  intractable 
problem.  Perhaps  as  useful  as  the  ability  to  perform 
such  analyses,  however,  is  the  extent  to  which  the 
armor  designer's  intuition  of  the  mechanics  of  penetra 
tion  is  enhanced  by  this  tool.  These  numerical  codes 
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are  easily  implemented  on  any  modern  computer  system, 
and  run  very  economically.  They  are  also  extremely 
amenable  to  user  modification  in  order  to  permit  easy 
implementation  of  various  constitutive  laws,  fracture 
models,  etc. 

Certain  areas  still  exist,  however,  for  significant 
improvement  in  this  treatment.  First,  the  codes  are 
presently  limited  to  zero-obliquity  impact  by  a 
projectile  whose  lateral  dimensions  are  small  compared 
with  the  region  of  influence  during  impact.  Relatively 
minor  code  modifications  would  be  necessary  to  include 
oblique  impacts  by  large  and  arbitrarily  shaped 
projectiles.  In  this  manner  the  influence  of  projectile 
geometry  could  be  modeled.  Second,  and  more  important, 
is  the  necessity  to  incorporate  more  acccurate  models 
of  material  response.  As  was  demonstrated  within  this 
report,  rather  sophisticated  constitutive  and 
fracture  algorithms  can  be  implemented  within  the 
codes  with  no  serious  difficulties.  More  work  is 
needed,  however,  to  determine  the  extent  to  which  these 
or  other  models  are  applicable  to  fabric  response  in 
the  ballistic  time  frame,  and  to  determine  the  numerical 
parameters  to  be  used  in  the  models. 

Two  examples  of  this  latter  problem  may  be  repeated 
here.  First,  recall  that  the  treatment  of  nonlinear 
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viscoelastic  effects  was  limited  not  by  the  code 
capability,  but  by  the  present  lack  of  understanding 
as  to  which  of  several  possible  nonlinear  constitutive 
laws  would  most  accurately  model  ballistic  response. 
Regardless  of  the  model  selected,  work  is  required  to 
obtain  experimentally  the  numerical  data  required  as 
input  parameters. 

Second,  the  otherwise  very  successful  materials 
parametric  study  described  in  the  report  is  deficient 
in  that  it  predicts  that  Kevlar  29  should  far  outperform 
Kevlar  49.  In  fact,  the  two  aramid  fabrics  perform 
almost  equally  well,  and  some  evidence  suggests  that 
Kevlar  49  is  actually  superior.  The  authors  have  no 

doubt  that  this  discrepancy  lies  not  in  the  wave- 
propagation  aspects  of  the  code  predictions,  but  in  the 
dynamic  failure  criterion  used.  Micrographs  of 
ballistically-fractured  fibers  show  extensive 
fibrillation,  and  evidence  exists  which  suggest  that 
Kevlar  29  and  49  differ  primarily  in  their  extent  of 
fibrillation  during  fracture.  Experimental  work  aimed 
at  elucidating  the  nature  of  the  fracture  mechanism 
is  needed;  incorporation  of  the  resulting  information 
into  the  penetration  codes  should  follow  without 
difficulty. 
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APPENDIX  A  -  The  FABRIC  Code 


General.  FABRIC  is  a  FORTRAN  coding  of  the 
numerical  analysis  of  fabric  impact  which  was 
described  in  Chapter  II.  In  its  present  .form,  the 
code  is  restricted  to  zero-obliquity  impact  on  an 
orthogonal  fabric  mesh  consisting  of  only  one  fiber 
type.  These  constraints  could  be  relaxed  through 
suitable  code  modifications.  The  code  was  developed 
and  implemented  on  MIT's  IBM  370/168  computer  system, 
but  was  later  implemented  without  difficulty  on  the 
NARADCOM  computer.  The  code  was  run  at  MIT  in  a  batch 
mode,  but  could  easily  be  modified  for  interactive 
terminal  operation:  this  would  likely  consist  primarily 
of  adding  terminal  queuing  for  data  input  and  graphical 
display  for  output  data. 

Code  input  and  output .  The  input  data  needed  by 
FABRIC  is  detailed  in  a  series  of  comment  lines  at  the 
beginning  of  the  code  listing.  Briefly,  these  include 
specification  of  the  impact  velocity  and  projectile 
mass,  the  fabric  idealization  (principally  the  number 
of  fibers  per  unit  length) ,  the  constitutive  and  frac¬ 
ture  properties  of  the  fibers,  and  such  run  paramaters 
as  maximum  alloted  time  and  printing  increment. 
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A  typical  data  input  set,  for  a  300-m/sec  impact 


on  a 

nylon  panel  is  given 

below: 

LINEAR  VISCOELASTIC 

FABRIC 

(BALLISTIC  NYLON) 

1.1 

300. 

19.533 

10.16  5280 

0.14 

0.0 

2. 

1.4142  1. 

3 

100. 

0.2 

5. 

36 

1  1 

Code  output  consists  first  of  a  series  of  values 
relative  to  the  initial  conditions  which  were  read  in 
and  which  the  computer  requires  in  order  to  begin  the 
recursive  calculations. 

After  each  time  increment  (or  less  often,  depending 
on  the  value  used  for  the  print  skip  increment  INC) ,  the 
code  prints  the  current  values  of  the  field  variables 
at  each  node  in  the  fabric  octant.  Currently,  these 
are  simply  dumped  in  order  of  the  calculation  scheme  as 
defined  by  Figure  10.  This  presentation  of  data  has 
been  sufficient  for  the  research  studies  discussed  in 
this  report,  but  for  production  design  work,  graphical 
or  some  other  high-order  output  would  likely  be  prefer- 
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EXECUTION  BEGINS ♦  ♦  ♦ 

1LINEAR  VISCOELASTIC  FABRIC  (BALLISTIC  NYLON) 


INPUT  PARAMETERS i 

YARN  DENIER ,  BENYRN  (DEN)®  0.52800E+04 

INITIAL  PROJECTILE  VELOCITY ,VPROJ  (M/SEC)®  0.30000E+03 

PROJECTILE  MASS*  PMASS  (GM)®  0 ♦  1 lOOOE+Ol 

FABRIC  PANEL  LENGTH  XL  (CM)®  O*  10160E+02 

ELEMENT  LENGTH  DXL,  (CM)®  0.29029E+00 

MAXIMUM  IMPACT  DURATION  TMAX  ( MICROSEC )®  0,20000E+01 

STABILITY  COEFFICIENT ,  CDTM®  0,14142E+01 

INCREMENTAL  TIME  DTM  < MICROSEC )=  0*69116E»-00 

NUMBERS  OF  LAYERS,  CLYR®  0 . 10000E+01 


STRAIN  WAVE  VELOCITY  CWAVE(M/SEC)=  0*29698E+04 
INITIAL  MODULUS  OF  YARN?  EYRN(GR/DEN>®  0 . 10000EL03 
BACK  UP  ELASTIC  SPRING  CONSTANT  XK  (GR/CM/CK)®  0*0 


NUMBER  OF  NODES  ALONG  MODEL  PANEL t  JT®  36 
NUMBER  OF  TIME  INCREMENT  STEPS ,  NTINC®  2 
PRINTING  FREQUENCY r INC-  1 


ACTUAL  FABRIC  MASS,FMASSA  (GH)  =  0*19533E+02 

MODEL  FABRIC  MASS,  FMASSM  (GM)  =  0»16689E+02 

CRIMP=FMASSA/FMASSM®  0  *  1 1704E+01 

UNIT  ELEMENT  MASS  *9* EOS  (UNITM)®  0,35877E+04 

HALF  UNIT  ELEMENT  MASS  #9.E  05  (HUNITM)-  0*17938E+04 


INITIAL  PROJECTILE  KINETIC  ENERGY,  XKE  (JOULE/GM)®  0.25342E+01 

MATERIALS  PROPERTIES ♦ 

OPTION-MATERIAL  MODEL,  IPT®  3 

INITIAL  MODULUS,  EYRN ( GR/DEN ) ®  0 ♦ 10000E+03 


VISCOELASTIC  MODEL. - STANDARD  LINEAR  SOLID  PARAMETER- 

MODEL  GLASSY  MODULUS  ( GPD ) ®  0* 10000E+03 

MODEL  VISCOUS  FRACTION®  0.20000E+00 

MODEL  RELAXATION  TIME (MICROSECONDS)®  0.50000E+01 


RUPTURE  STRAIN  OF  INDIVIDUAL  FIBER, EPSB®  0.14000E+00 
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Program  Requirements.  In  its  present  form,  the 
FABRIC  code  is  self-contained,  needing  no  external 
software  support.  All  its  subroutines,  for  instance, 
are  contained  within  the  listing  given  below.  One 
qualification  to  this  statement  however,  is  the  inclu¬ 
sion  of  a  call  to  Subroutine  UERTST  within  the  nonlinear 
equation-solving  subroutine  ZNOLNR.  UERTST  is  an  error¬ 
handling  routine  available  through  the  proprietary 
"International  Mathematics  and  Statistics  Library" 

(IMSL) .  If  IMSL  is  not  available  at  the  user's  location, 
the  call  to  UERTST  could  be  removed  with  little  risk: 
no  error-handling  capability  was  needed  during  any  of 
*  the  computer  experiments  described  in  the  body  of  this 
report. 

FABRIC  needs  no  tape  or  disk  support;  all  compu¬ 
tations  are  performed  directly  within  core.  Core 
requirements  and  run  times  are  dependent  both  on  problem 
specifications  and  on  the  computer  system,  but  by  way 
of  illustration  some  parameters  observed  during  a  typi¬ 
cal  run  on  the  IBM  370/168  system  will  be  mentioned. 

During  an  impact  simulation  of  a  Kevlar  29  single  layer 
at  300  m/sec,  the  computed  time-to-penetration  was 
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42.87 


sec.  The  panel  was  idealized  as  having  a 
fiber  spacing  of  0.3175  cm,  and  the  time  step  for 
optimum  code  stability  was  0.32234  ^Csec.  A  total  of 
133  time  increments  was  therefore  necessary  for  the 
run,  with  each  time  increment  involving  an  additional 
node  relative  to  the  previous  step.  The  total  run 
time  for  this  job  was  1.089  minutes  ($12.92  at  weekend 
rates) ,  and  a  total  of  182  kilobytes  of  CPU  core  was 
required.  When  this  same  problem  was  run  at  NARADCOM, 

3  min  50  sec  of  CPU  time  was  required,  which  illustrates 
the  system  dependency  of  such  job  parameters. 


FLOW  DIAGRAM  AND  LISTING 


lift 
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Initialize  Arrays: 


compute  new  velocities 
Eq.  20-21 

V 

compute  new  coordinates 
Eq.  '? 

1 

compute  element  strain 
Eq.  23-24 

3 

compute  element  tension 
(subroutine  TENSN) 


Eq.  25,40,46 


compute  kinetic  and 
strain  energies 
_ Eq.  29,30 
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FILE:  FABRIC  FORTRAN  A  CONVERSATIONAL  MONITOR  SYSTEM 
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APPENDIX  B  -  The  XOVER  Code 


General .  XOVER  is  the  FORTRAN  program  written  to 
carry  out  the  calculations  described  in  Chapter  IV 
of  this  report.  It  is  identical  in  concept  to  the 
FABRIC  code,  differing  primarily  in  the  dimensionality 
of  the  problem.  XOVER  is  essentially  a  treatmei  >.  of 
transverse  ballistic  impact  on  a  single  fiber  as 
described  earlier  by  Roylance  [18] ,  but  it  has  been 
modified  to  include  a  second  fiber,  transverse  to  the 
first,  which  receives  its  loading  from  the  motion  of 
the  first  (primary)  fiber.  The  code  is  made  more 
complicated  than  the  single-fiber  case  by  the  necessity 
of  allowing  motion  in  all  three  directions  for  the 
secondary  fiber,  and  in  computing  the  initial  values  of 
the  secondary  fiber  in  terms  of  the  primary  fiber 
motion.  As  described  in  the  text,  allowance  is  made 
for  slippage  of  the  secondary  fiber  along  the  primary. 

Code  requirements.  As  was  true  for  the  FABRIC 
code,  XOVER  requires  little  hardware  or  software 
support.  At  present,  the  code  includes  a  call  to  an 
MIT  library  subroutine  PRTPLT  for  the  purpose  of 
obtaining  a  rough  plot  of  the  nodal  variable  values 
on  the  system  line  printer.  This  subroutine  call  could 
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be  removed  without  affecting  the  performance  of  other 
parts  of  the  code.  For  a  typical  run,  in  which  two 
crossed  Kevlar  fibers  were  impacted  at  300  m/sec,  the 
job  run  time  was  0.305  minutes  and  170  kilobytes  of 
core  was  reserved. 

Definition  of  program  variables.  XOVER  is 
conveniently  described  by  means  of  a  listing  of  the 
principal  program  variables:  (*  denotes  input  data) 


A 

CROSS* 

CZ 

DEN* 

DL 

DT 


8.826  x  10°  (dt/dl) 

Initial  crossover  position  (fraction  of  XLl) 
Wavespeed  (m/ sec ) 

Fiber  denier 
Length  increment  (cm) 

Time  increment  (sec) 


El ( J) ,  E2(J)  Strain  at  jth  node  in  primary  and  secondary 
fiber 


EIOLD(J) , 
E20LD(0) 

EMAX* 

EEMAX 

G* 

IPLOT* 

11,12,13* 


Strain  at  previous  time  increment 

Maximum  strain  permitted 

Maximum  strain  in  fibers 

Instantaneous  modulus  (gm/den) 

.EQ.  0  if  no  printer  plot  is  desired 

Node  numbers  at  which  plot  desired  (13  on 
secondary) 
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JX 

JJX 

KMAX 

LPSKIP* 

LTSKIP* 

NLINC* 

PMASS* 

PVELOC* 

SLIDE* 

TAU* 

TITLE* 

T1(J)  , 

T2  (J) 

TIOLD(J) , 
T2)LD(J) 

TMAX* 

U1(J)  , 

U2  ( J) 

VI  (J)  , 
V2(J) 


Crossover  node  on  primary  fiber 
Initial  crossover  node 
Maximum  number  of  time  increments 
Length  print  skip 
Time  print  skip 

Number  of  length  increments  on  primary 
fiber 

Projectile  mass  (gm) 

Projectile  velocity  (m/sec) 

Slip  factor  (1-no  slip,  0-no  friction) 
Viscous  relaxation  time  (sec) 
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PMASS , 
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2E  10. 
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IPL0T, 

11,  12,  13 

Typical  Output.  Typical  output  from  the  XOVER  code 
for  impact  on  two  crossed  Kevlar  29  fibers  at  300 
m/sec  is  given  below: 
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